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Abstract 

The Casimir repulsion between a metal and a dielectric suspended in a liquid has been thoroughly 
studied in recent experiments. In the present paper we consider surface modes in three layered 
systems modeled by dielectric functions guaranteeing repulsion. It is shown that surface modes 
play a decisive role in this phenomenon at short separations. For a toy plasma model we find the 
contribution of the surface modes at all distances. 
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I. INTRODUCTION 



The existence of repulsive Casimir forces follows straightforwardly from the Lifshitz the- 
ory [l| for a special choice of the materials. The phenomenon has been discussed in several 
theoretical papers mainly for materials with nonunit magnetic permeability. 

Recent advances in the design of metamaterials 8| demonstrating nontrivial magnetic per- 
meability at optical frequencies stimulated experiments seeking after the Casimir repulsion. 
Though some experimental groups have reported noticeable decrease of the attraction 9| in 
the presence of metamaterials, to our knowledge, nobody has succeeded so far in reversing 
the sign of the force. 

The repulsive Casimir forces between ]5urely dielectric materials now appear to be more 
promising for observation and might turn out useful for nanomechanical systems. There 
repulsion is achieved by filling the space between the bodies with a medium with a precisely 
chosen dielectric permittivity. 

It was experimentally shown [lO, 
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121] for short separations (<50 nm) that repulsive 



van der Waals forces between specially chosen dielectrics with intervening liquids of low- 
polarity (eg. cyclohexane , ethanol, or bromobenzene) agree with theoretical predictions of 
the Lifshitz theory including the retardation contribution. The repulsive Casimir forces at 
large separations u p t o several hundred nanometers have not been experimentally studied 
until very recently 

m Q, Q, Q, Q 



The setup is the following. Two surfaces '1' and '3' with dielectric permittivities ei and 
£3 are separated by a gap filled by medium '2' with dielectric permittivity £2- In recen; 
experiments or proposals medium '2' is a liquid, ethanol 16, 17] or bromobenzene |l8| 
while the two surfaces are made of silica and gold. In order to obtain a repulsive force the 
respective dielectric functions have to satisfy the relation 



(1) 



in the frequency range relevant for the force measurement. 
The Casimir force is given by 
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Here L is the width of the gap between the plates; rf2(ico', k), p = TE, TM, are the reflection 
coefficients at imaginary frequencies for the surfaces facing the medium '2' 

r™M = ^^^^i^, rf/M = -^^, (3) 

e2Ki + eiK2 Ki + K2 



with Kj{iuj) = ^Jk?' + e j{iuj)uj'^ / , j = 1, 2, 3. 

In the present system the Casimir force is repulsive. Indeed, at short distances, L « Xp^i, 
the force is approximated by F ^ —Hi23/3L^ . Its magnitude and sign is defined by the non- 
retarded Hamaker constant [lol 

oo 

= ^ I du^l^ . (4) 

with Aj2 = {£2ii^) — gjfici; 



57r^ I n-^ 
"=i 



constants are negative 
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)/{e2{iuj) + ej{iuj)). Under the condition ([T]) the Hamaker 
16| and the force is repulsive. 
At long distances L » Xp i, where the largest contribution to the force comes from small 
frequencies and small wave vectors k, the reflection coefficients tend to their the static values 



rJ/.r™«l„„«Z-^^, , = 1.3. (5) 



Plugging 



tions 



them into (Ej) one gets a rough estimation of the force at large plates' separa- 



20l |. It is easy to check numerically that if the materials obey ([T]), especially at low 



frequencies, the force is repulsive in the long distance limit. These considerations help to 
choose the materials for the devices based on repulsive Casimir forces Q, [3]. However the 
origin of this repulsion itself is not quite clear. 

In the following we will consider this question in relation with surface modes. The 
electromagnetic quantum fluctuations which give rise to the Casimir effect obey the Maxwell 
equations. The corresponding boundary-value problem has two types of solutions: the 
propagative waves, which will be called photonic modes in the following, and the waves living 
on the interfaces and exponentially decaying outwards. These surface modes exist only in the 
TM polarization. Further on we call them surface plasmons or polaritons depending on the 
model which describes the material. The term "plasmon" is reserved for the plasma model. 
The expression ([2]) comprises the contributions from both propagative and surface modes. 
The present paper pursues their respective roles in the Casimir repulsion by considering a 
model system consisting of three layers of dielectric material the plasma frequencies of which 



are chosen such that condition [T] is fulfilled. Section II and III describe the definition of the 
surface modes and their contribution to the Casimir energy. In section IV we perform a 
calculation of the repulsive Casimir force for the recent experiment ISj where we use as 
dielectric functions the two oscillator model for the silica surface and the bromobenzene 
liquid filling the gap. The second surface is covered with gold for which we use the Drude 
model as dielectric response. A comparison between this calculation and the predictions of 
our model system allows to define the limitations of the latter. The paper finishes with some 
conclusive remarks in section V. 

We will assume that the dielectric properties of the liquid do not change as its layer 
becomes thicker or thinner. We will not consider the hydrodynamics of this system either. 
Of cause, a body moving in a liquid due to the Casimir attraction or repulsion experiences 
drug force, which depends on the velocity. This hydrodynamic force should be taken into 
account in the experiments. Luckily it does not depend on the reason of the movement, and 
may be estimated separately jl4, Q, [l8 |. 



II. THE INTERACTION OF THE SURFACE PLASMONS 

Let us first formulate the interaction of the surface modes to understand the nature 
of Casimir repulsion. For simplicity we consider materials described by plasma model, 
Ei = 1 — ujpju^, i = 1..3, where Upi is the material's plasma frequency. The frequency of 
the single surface plasmon living on the interface of medium i, i = 1, 3, with medium '2' is 
given by 



sp 
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^ 2k\' + + - Mk^c^ + {u^l, - ulf ^ (6) 



When ^ the single surface plasmon frequency tends to ujp2, provided ujpi/ujp2 < 1, 
otherwise it approaches Upi. In the limit /c — > oo it tends to y^Ci^^~+~ci^/v^, i = 1, 3. 

The surface modes evanescent in the direction of the gap between the plates are coupled 
through the cavity. Their frequencies are defined by the equation 

with Qi = \/k'^ — EiLu'^/c^. Let us compare the behavior of the solutions of Eq. ([7]) in two cases 
which result in opposite signs of the Casimir force. Fig.l shows the plasmon modes inside 
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FIG. 1: The plasmon modes in two set-ups: (a) ujp2 < ujpi,ujp-^ (attractive Casimir force), (b) 
Ljpi < ujp2 < ujp3 (repulsive Casimir force). The sohd hnes correspond to coupled symmetric and 
antisymmetric plasmons. The dashed lines show the single surface plasmons living on the interfaces 
1-2 and 2-3. The dotted lines mark the boundaries between propagative and evanescent sectors. 

the "sandwich" (a) Ljp2 < Upi,ujp3 (attractive Casimir force), (b) Upi < ujp2 < ujp3 (repulsive 
Casimir force). The coupled plasmon modes are plotted as solid lines. The dashed hnes 
correspond to the plasmons living on the single interfaces. Situation (a) corresponds to 
well known case where the gap in-between the two plates is filled with vacuum fluctations 



{ujp2 = 0). It has been studied in detail in the papers [22|, |23|. For the sake of comparison 
with case (b) let us briefly recall the results. Fig. la shows the solutions of equation ([7]) when 
a = Up2/u!pi = 0.5, P = Ups/ujpi = 1 (equal slabs interlaid with material having ujp^2)- The 
dotted lines starting from u/upi = a and /3 separate the propagative and evanescent sectors 
respectively in the gap and in the slabs. We denote the solutions of ([7]) hj u!±, where is 
usually referred to as antisymmetric plasmon and u;_ as symmetric plasmon. The latter one 
lies entirely in the evanescent sector , k"^ —u"^ /c^ + 00^2/0"^ > 0. In contrast, the antisymmetric 
plasmon penetrates into the propagative sector if k < 



P+ = — \ \ ce } , (8) 

where A = Lupi/c is the dimensionless distance. The coupled plasmons u!± surround the 
single surface plasmon solution u^^, and u;+|fe_»o = ^pi, ^-\k->Q = ^p2- At large wave vectors 
u;Jf, uj^2- For equal slabs, meaning /3 = 1, these limits coincide (dashed line). 
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The modes which are evanescent in the gap are evanescent within the slabs as well, because 
on Fig. la the borderline between the evanescent and the propagative sectors in the gap lies 
below the corresponding borderline for the slabs. 

Fig. lb shows just the opposite situation. The dotted lines starting from uj/ujpi = 1 and 
a separate the propagative and evanescent sectors respectively in the gap and in the slab 
with ujpi. Now the borderline between the two sectors referring to the gap lies above the 
one corresponding to the slab. In the area between these curves the equation ([7]) has no real 
solutions as g|, qj > 0, but qf < 0. Solving ([7]) with k"^ — u'^/c^ + 1^^2/0^ 0, one finds that 
the mode a;+ exists for k > k^^ 



k^ = ^ 



1 + / 



(9) 



^ V a^ 1 ^^^^ ^ 

The coupled plasmons uj± lie inside the area enveloped by the single surface plasmon solutions 
a;32 and defined by ([6]), and uj-\k^o = Upi, uj+\k^Q = Up2. 



III. THE VACUUM ENERGY OF THE SURFACE MODES 

In the following we will concentrate on the Casimir energy contribution of surface modes 
in the case of a repulsive force. The renormalized vacuum energy of the interacting surface 
plasmons living on the plane mirrors is formally given by 

00 

^"-|e/^i-jl.. (10) 

where lim/^_^oot^+ = ^12; li^iL^oo'^- = "^32? ^- = 0; fc+ is defined by ([9]). 

To calculate the surface plasmon energy, we first introduce the dimensionless variables 
K = kL,n = ujL/c, and ^Ipi = A, = aA, ^Ip^ = f3A in (iTllTOll 

{CO 00 \ 

J dKKisi^ - nil) + J dKKin. - nZ) \ . (n) 
K+ J 

Then we change the integration variable in ffTTl) . K ^ Q = ^K"^ — n? and write the 
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FIG. 2: The normalized vacuum energy: (a) symmetric plasmon lo^ (dashed hne), antisymmetric 
plasmon a;+ (dash-dotted hne), sum of both (sohd hne); (b) plasmon energy (dashed line), photon 
energy (dash-dotted line), total Casimir energy (solid line). ujp2/ujpi = SjUps/wpi = 25. 

renormalized energies of the symmetric and antisymmetric plasmons as 

{oo 

-A2 ^ 

-II dQ'nil + l{nl-nl,)\^\, (12) 

where X = Kl-nl^{K+). 

The numerical results for the plasmon energy are given in FigEJ^a). The distance is nor- 
malized by the largest plasma wavelength in the system, Api = l-ncjiOpx. The plot shows the 
reduction factor which is the ratio of the plasmon vacuum energy and the Casimir energy 
of perfectly conducting plates, Ec = —hcir^ /720L^. As the Casimir force between perfectly 
conducting plates is attractive, a negative or positive reduction factor corresponds to repul- 
sion or attraction respectively. For the energies, the arguments goes in the opposite way. 
The energy of the antisymmetric plasmon a;+ is positive, corresponding to repulsion with a 
negative reduction factor, while the energy of the symmetric plasmon U- is negative yielding 
attraction and thus a positive reduction factor. The antisymmetric plasmon dominates at 
short distances, while the symmetric plasmon makes decisive contribution at larger separa- 
tions. Thus the total plasmon interaction is repulsive at short separations and attractive at 
medium and long distances. Here we estimate the length scale with respect to Api. 
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In [22|, |23|] , it was shown that for two plates separated by vacuum the attractive Casimir 
energy is a result of cancelations between plasmon and photon contributions at all distances. 
Moreover the plasmon and photon energies have different signs and decrease more slowly 
with distance (~ L~^^^) than the total energy (~ L~^). 

Fig|2](b) shows the reduction factor for the total Casimir energy when condition ([T]) is 
met. Clearly it is negative at all distances corresponding to a repulsive Casimir interaction. 
It comes out as the sum of negative photon contribution and the plasmon contribution, 
plotted in Figl2](a). The plasmon contribution dominates at short separations. At medium 
and long distances the photon contribution prevails assuring repulsion. At distances large 
with respect to Api the energy decreases as ~ Ei(l,2Q;A), A = 27rL/Api. The decrease is 
much steeper than in any set-up yielding attraction. 

The attraction of the surface plasmons at short ranges explains the sign and the mag- 
nitude of the Casimir force between plates separated by vacuum (or any material with 
u!p2 < ujpijUJps). In the following we will show that under the condition ([T]) the surface 
plasmons produce a repulsive contribution at short distances. 

The equation for the interacting surface plasmons may be solved explicitly for large wave 
vectors which correspond to short distances: 

1 



,2 \2 



Then the vacuum energy is 



where 



+4(^=,-^J,,)(^=3-^JJe-=l'l'-l'/'}. (13) 





^± = {2«2 + l + /32^ [(1-/32)2 

If at large frequencies condition ([1]) is satisfied, then > yielding repulsion. The simplest 
particular case is /5 = q;2 > 1 {ijjp^2/^p,i = ^p,2,/^p,i) corresponding to equal ratios between 
the plasma frequencies. For large a, ip{a oo) 0.67q;. If a = 1-1-5, Y{5 ^ 0) ~ 5^/2. For 
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example, 1^(1.1) = 0.00497. Thus, the vacuum energy of the interacting surface plasmons is 
positive corresponding to repulsion. 



IV. REPULSIVE CASIMIR FORCE USING THE DIELECTRIC RESPONSE OF 
THE MATERIALS 

Let us now compare the repulsive force due to the interaction of the surface modes ob- 
tained in the plasma model with the corresponding total Casimir force calculated according 
to the zero temperature Lifshitz formula ([2]). On Figl3]we plot the force ([2]) normalized by the 
Casimir force between perfect conductors as a function of the dimensionless distance L/Api, 
Api = 27rc/co'p_i. The calculation is done for tree layers with = Woj^^juP' , a = ujp2/ujpi = 5, 
f3 = ujp^/ujpi = 25 (solid line). At short distances the force coincides with the short distance 
asymptote of the interacting surface plasmons, F^f = —dE^^/dL, F^^/Fc = — 7.38L/Api 
(dotted line). At long distances (L >> Ap2) the force is repulsive and decays much faster 
than the force between two plates separated by vacuum, F ~ exp(— 2co'p2-L/c)/L, where u!p2 
is the plasma frequency of the material filling the gap. This behavior is explained by the di- 
vergence of £2 for C(j — * within the plasma model. If dielectric function of the intermediate 
material '2' is finite at low frequencies, the reduction factor saturates at large distances. 

An example of a more realistic set-up is the system Silica-Bromobenzene-Gold. Repulsive 
dispersion forces in a similar system were measured by {l^ at short distances, the results 
appeared to be consistent with the calculated non-retarded Hamaker constants. Precisely 



the same system was studied in a recent experiment [18[ at separations from 20 nm to 
several hundred nanometers. Here we do not consider separations larger than 300 nm, so 
that we need not to account for temperature corrections. In the following we first calculate 
the reduction factor applying the zero temperature Lifshitz formula ([2]) with the materials 
described by Drude and Lorentz models. 

For gold eAu{i^) = 1 + ^Au/l^i^ + ^)]> where uau = 1-367 • 10"^^rad/s, 7 = 5.316 • 
10^'^rad/s [2]. To simplify the analysis for bromobenzene and silica we confine ourselves to 
two oscillators in the standard multiple oscillator model 

^^(-) = 1 + ,^^7^ ^2 + ,^^7^ (16) 

1 + (^M^)2 1 + {uj/uj}jyy 

where i = 1,2 and correspond to silica and bromobenzene respectively. The various param- 
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eters are C)^ = 0.829, 00} = 0.867 ■ lO^^rad/s, C^y = 1.098, ooljy = 2.034 ■ 10^^ rad/s [l9|; 
C]j^ = 2.967, = 5.47 ■ 10%ad/s, Cfjy = 1.335, to^y = 1.286 • 10^^ rad/s [10]. With 
these values of the parameters the condition ([T]) is satisfied for < 9 • lO^^rad/s, where 
£si02 < ^CaHsBr < £au- For > 3 ■ 10^*^rad/s, Eau < ^CeH.Br < £si02- Both regions 
contribute to repulsion within the two oscillator model for bromobenzene and silica. The 



applicability of this model is discussed in [l7|, [l8|. 



In the present situation the solutions of ([7]) are interacting surface polaritons. At close 
separations realized by a thin layer of bromobenzene, where the decisive contribution to the 
force comes from large wave vectors and high frequencies, the absorbtion is negligible, and 
the models used for the materials are reduced to the plasma model with effective plasma 
frequency Wpf^ = [C^yW^y + »^ = 1^2, ujps = uau- We find w^g^ = 2.131 • 

lO^^rad/c, uj^'J^^HsBt = 1-488 ■ lO^^rad/c. 

At close separation the surface polaritons turn into surface plasmons, and the for- 
mulae (113tll5p for the the interacting surface plasmons become valid. But one should 
keep in mind that they describe correctly only short distance regime. Substituting a = 
^p"/^p" = 0.698, /3 = ujps/uf/ = 0.641, we get F^p/Fc = -0.03355 L/A;{^, where 
Apf'^ = A^'^Q^ = 27ic/ujp{-^ = 88.44nm is the effective plasma wavelength for the composed 
system. In this regime the reduction factor for the interacting surface modes is negative 
yielding repulsion. It coincides with the total reduction factor for the force at short dis- 
tances. 

At medium and low frequencies the dielectric functions of the materials are complex, 
and the calculation of the vacuum energy corresponding to the surface modes is not so 
straightforward. The total reduction factor for the force is shown as the dashed line in 
Fig. 3. Clearly repulsion sets in for distances of the order of or larger than the effective plasma 
frequency while the toy model gives repulsion essentially for distances much smaller than that 
and gives a negligible force for larger distances. In contrast at large distances our calculation 
using the two oscillator model yields a negative force with a reduction factor which reaches 
T] ^ —0.057. Within the plasma model the reduction factor decays exponentially, while 
within the two oscillator model it tends to a constant. This shows that the repulsion at 
large distances does not have its origin in the interaction between surface modes but stems 
from the repulsive contribution of propagating modes. 

The dotted lines are the short distance plasmon approximations to the model system 
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FIG. 3: Reduction factor rjp = F/Fc at short plate separations as a function of dimensionless 
distance A = L/Xpi; a = Wp,3/wp,2 = ^p,2/^p,i- The dotted hues are the respective short distance 
asymptotes of the surface plasmon interaction, a = ujp,3/iOp^2 = ^p,2/^p,i- For the system Si02- 
CeHsBr-Au Api = iT^cju^l^l^Q^ = 88.44nm. 

and for the realistic dielectric response functions. For the latter the plasmon approximation 
turns out to be valid for distances up to twice the effective plasma wavelength, that is about 
170 nm. Indeed, the present calculation does not pretend to be a precise description of the 



recent experiment [18| , but the system Au — C^H^Br — Si02 is used as an illustration. 



The calculation of the Casimir force at zero temperature for the solid-liquid-solid system 



using measured dielectric functions of all involved materia 
millimeters down to subnanometers was carried out in 
neously to the present paper. 



s for the wavelength range from 



24| . which was submitted simulta- 



V. CONCLUSION 



In the present paper we have studied the Casimir force in the systems which meet the 
condition ([1]). The magnitude of the force in the present setup is small. For example, in the 
system 5*^02 — CqH^Bt — Au, with the 150 nm layer of bromobenzene, it is only about 5% 
of the force between gold covered plates in vacuum. 



From previous studies 



22 



231] we remember that when two plates in vacuum form a 



cavity, the total energy of the coupled surface modes is negative, corresponding to attraction. 
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We have shown that within the plasma model the surface modes repel at short distances when 
the materials satisfy the condition ([1]). Moreover the Casimir repulsion is then completely 
due to the repulsion of the surface modes. At medium and long distances the interaction 
of the surface modes becomes attractive, but the dominating repulsive contribution of the 
propagative modes leads to a total repulsive force. 
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